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ABSTRACT 

Exact bolometric light curves of supernova sho ck breakouts are derived based on the universal, non relativis- 
tic, planar breakout solutions ( Sapir et al.ll201 lab , assuming spherical symmetry, constant Thomson scattering 
opacity, k, and angular intensity corresponding to the steady state planar limit. These approximations are ac- 
curate for progenitors with a scale height much smaller than the radius. The light curves are insensitive to the 
density profile and are determined by the progenitor radius R, and the breakout velocity and density, Vq and 
po respectively, and k. The total breakout energy, Ebo, and the maximal ejecta velocity, v,nax, are shown to be 
Ebo = 8.07r7?^K"'cV() and Vmax = 2.0vo respectively, to an accuracy of about 10%. The calculated light curves 
are valid up to the time of transition to spherical expansion, fjph « /?/4v(). Approximate analytic expressions for 
the light curves are provided for breakouts in which the shock crossing time at breakout, to = c/ npov^, is^R/c 
(valid for R < 10^* cm). Modifications of the flux angular intensity distribution and differences in shock arrival 
times to the surface, Afasym, due to moderately asymmetric explosions, affect the early light curve but do not 
affect Vmax and Ebo- For 4vo ^ c, valid for large (RSG) progenitors, L oc f""^/^ at max(Afasym,^/c) <t< fsph 
and R may be accurately estimated from R^2x 10'-'(L/10'*^ ergsec"')^''^(f/lhr)^/'^. 
Subject headings: radiation mechanisms: non-thermal, shock waves, supernovae: general. X-rays: bursts 



expeciea lo oe e r mueg (n^oigatei iiv/^ iraua 
& Chevalier 119781: lEnsman & BurrowTll992l: 
[cKea,1999: ^Blinnikov et alj 12000: iKatz et alj 



1. INTRODUCTION 

During a core collapse supernova (SN) explosion, a strong 
radiation mediated shock (RMS) traverses the exploding 
stars' mantle/envelope. Once the shock reaches the sur- 
face of the star, a burst of high en ergy radiation (UV 
to 7-rays) is expected to be emitted ([Colgatel 119741: iFalkl 
[1978: Klein, 
Matzner & McKee 

11009; Piro et al.iSoiOtlMakar & SariJ|201( 

The observed properties of this breakout were derived in 
earUer analyses using either anal ytic order of magnitude 
estimates (e.g. Matzner & McKeeJpH iKatz et alj 120091 : 
iPiro et all 120 lOt .Nakar & Sari. 2010b or numerical calcu- 
lation s for particular progenitors (e. g. Ensman & Burrows 
1992'; 'Blinnikov et al.' ^OOff; 'Utrobin^007; Tominaga et al. 
2009; Tolstov 2010; Dessart & Hillier 2010; Kasen et al. 
201 II) . Here we provide an accurate description of the time 



dependent radiation emission following shock breakout for a 
general progenitor without an opti cally thick wind. 

In the first paper of this series (S apir et aDl201 la l. hereafter 
Paper I) we have solved the problem of a non-steady planar 
RMS breaking out from a surface with a power-law density 
profile, p (X x", in the approximation of diffusion with con- 
stant opacity (a similar solution for exponential atmospheres, 
va lid only for the early p art of the planar breakout, was given 
in iLasher & Chanlll979l) . In this paper we use the results of 
the planar calculation to derive the observed bolometric prop- 
erties of SN shock breakout bursts, taking into account limb 
darkening. In a third paper (ISapir et al.ll201 Ibl) we calculate 
the temperature profiles and spectral properties of the burst 
assuming local Compton equilibrium and photon generation 
by Bremsstrahlung (see § 5 of paper I). 
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The derivation of exact light curves for a general progeni- 
tor is possible thanks to the universality of the planar breakout 
solutions (Paper I). Radiation escapes the shock front, produc- 
ing the observed breakout burst, when the optical depth r of 
the plasma lying ahead of it is equal to the optical depth of the 
shock transition layer, Tsh = c/vsh (Reaver 1976). We denote 
the shock velocity and (pre-shock) density at this point by vo 
and Po respectively (see §|2]for exact definitions). Measur- 
ing length, time and mass in units of xq = c/npoVo, to = xq/vq 
and nio/xfi = c/kvq respectively, the RMS breakout solution 
is universal and depends only on the density power law in- 
dex n. It was numerically found that the luminosity depends 
weakly on the density structure. In fact, throughout most of 
the emission, the luminosity changes by less than 25% for n 
in the range 1-10. The fact that the planar luminosity curve 
changes so little for a wide range of density power law in- 
dexes suggests that the results are insensitive to the decreas- 
ing density structure, and are applicable to profiles which are 
not power laws. 

The results presented in this paper, based on the non- 
relativistic planar breakout solution, have corrections of or- 
der Vo / c. Estimating these and higher order corrections re- 
quire a relativistic calculation of a planar shock breakout, 
and are beyond the scope of this paper We do estimate one 
obvious first order correction resulting from the transforma- 
tion of the comoving calculated flux to the observer frame in 
§ |B] An additional complication occurs at shock velocities 
v/c ^, 0.3, where the post shoc k temperatures reach 50 keV 
( IWeaverl 119761: iKatz et all 12009 ) and electron-positron pairs 
are created and increase the opacity of the material. Our dis- 
cussion is applicable to smaller velocities/lower temperatures. 
Order of magnitude estimates for relativistic shock breakouts 
can be obtain ed by using the steady state solution of a rela- 
tivistic shock (iBudnik et al.ll2010tlKatz et alJl2009h . 

This paper is organized as follows. The planar non- 
relativistic RMS breakout solution is shortly described in § |2] 
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The application of the planar solution to a SN breakout is 
discussed in § [3] The planar breakout parameters vq and po 
are expressed in terms of progenitor properties and explosion 
energy, and complications due to spherical expansion, asym- 
metry and relativistic corrections are discussed. Numerically 
calibrated analytic expressions for the total breakout energy, 
Eqo, the asymptotic velocity of the fastest moving ejecta, 
Vmax, and the luminosity at late times, t > max(Afasym,^/c), 
are given in § |4] These are all insensitive to small deviations 
from spherical symmetry and from the constant flux angular 
intensity distribution. In § |5] the observed bolometric light 
curve during breakout is calculated, assuming that the shock 
arrives at the surface simultaneously and taking into account 
light travel time effects. Approximate analytic expressions for 
the light curve are derived. Some of the first order corrections 
in v/c to the burst properties are analyzed in § |B] Our results 
are compared to those of previous studies in § |6] The main 
results and conclusions are summarized in § |7] Appendices 
§Ia1-§Ic1 provide details omitted from the manuscript. 

2. PLANAR NON-RELATIVISTIC RMS BREAKOUT 

The problem of a non relativistic RMS breaking out from 
the surface of a plan ar decreasing density profile was solved 
in Paper I (see also iLasher & Chanlll979l) . The analysis is 
preformed by neglecting the thermal energy of the matter and 
by approximating the radiation energy transport by diffusion 
with constant opacity k. 

The initial density profile is assumed to be a power-law, p cx 
x", where x is the distance from the surface. The initial density 
and the asymptotic shock velocity at large optical depth are 
parameterized by 



«/(«+!) 



and 



p(T) = po(vor/c) 



v,(t) ^ vo(vor/c)-''«"/(«^" 



(1) 



(2) 



respectively, where t = kJ pdx is measured with respect to 
the surface and where /3„ is a function of n, which can be ob- 
tained num erically by so lving the pure hydrodynamic shock 
evolution (Sakurai 1960). The values of po and Vq are the den- 
sity and velocity at the point t = c/vsh that would have been 
obtained in a pure hydrodynamic shock propagation. 

The limits n — > and n — > oo are well defined and corre- 
spond to a homogeneous density distribution and an exponen- 
tial profile dHavesll 19681) . 



P- 



pq exp(K/cio Ax), Ax = x(p)-x(po) 



(3) 



respectively. The value of f3„ in Eq. ^ decreases monotoni- 
cally from 0.207 to 0.176 for < n < oo. 

As mentioned in § [T] measuring length, time and mass in 
units of xo = c/kpoVo, k) = xq/vq and mo/xg = c/kvo respec- 
tively, the RMS breakout solution is universal, i.e. depends on 
the dimensionless parameter n alone. In particular, the planar 
emitted energy flux, i.e. the luminosity per unit area, can be 
expressed as 

Cit) = Coc(^-^,ny (4) 

where Co = poVo and C = C/Cqis given in table 3 of Paper I 
for « in the range 1-10. The shock crossing time at breakout 
is defined as 

fo = — (5) 



The reference time fief corresponds to fpeak in Paper I, the time 
at which the luminosity per unit area C peaks in the planar 
solution. Since the observed SN light curves do not peak at the 
same time (due to light travel time effects) we have replaced 
fpeak with fref to avoid confusion. 

As can be seen in table 3 of Paper I, the luminosity depends 
weakly on the density structure. In fact, throughout most of 
the emission, the luminosity changes by less than 25% for n 
in the range 1-10. The fact that the planar luminosity curve 
changes so little for a wide range of density power law in- 
dexes suggests that the result is insensitive to the decreasing 
density structure, and likely represents also profiles which de- 
viate from a power law structure. 

For convenience, we include a surface area of 4ttR^ in the 
expressions below. The luminosity, L(f ) = AirR^Cit), may be 
written as 



Lit): 



to 



where Lq is the breakout luminosity defined by 

Lo- 



; 47r/?VoVo. 



(6) 



(7) 



At late times, t ^ tg, th e luminosity follows L(t) cx f"^/^ 
(iPiro et al.l2010tlNakar & Sari 2010fl InPaperl it was found 
that the exact solutions for the luminosity and its integral, 
E{t) = J L{t')dt' , are well approximated by 



L(f) = Lo 



E(t) = E, 



to 



t - fref 

a,to 



-i/i 



(8) 



where 



: 2.0 X 4ttR 



2V0C 



Loo = 0.33 x 47r/?VoVo, 

a, = (3Loofo/£oo)' = 0.125. (9) 

Eqs. ([8]l and (|9]l describe the emitted flux to an accuracy of 
better than 30% (10%) in L(f) (E(t)) for 1 < « < 10 and 1 < 
(f-fref)/fo < 100. Individual fits to different values of n allow 
higher accuracy. Note that £00 is the total energy emitted in 
the planar approximation. 

Finally, in the non relativistic approximation in planar 
geometry an exact relation exists between the velocity of 
the outermost mass element and the emitted luminosity 
(ILasher & Chanl[T979t iSapir et al.ll20TTah . 



v(f)=. 



at')dt' = 



KE(t) 

AttR-c' 



(10) 



In particular, the asymptotic value of the velocity of the sur- 
face is 

Equation ( fTOl l simply states that photons that hit a given par- 
ticle transfer all their momentum to the particle on average . It 

^ for a constant density profile, n = 0, the decline is steeper, L{t) oc 
ISapir etal.120 Hal) 
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holds for any elastic scattering which has forward/backward 
symmetry, regardless of whether the diffusion approximation 
is valid or not. 

3. APPLICATION OF THE PLANAR SOLUTION TO A 
SN BREAKOUT 

In this section the application of the solution of shock 
breakout in the planar ap proxi mation to supernova shock 
breakouts is discussed. In § 13. li the relation between the pro- 
genitor parameters and the breakout parameters is briefly re - 
viewed. The effects that need to be taken into account in 
calculating the l ight curves using the planar solutions are de- 
scribed in § 13.21 The limitations of the planar approximation 
are discussed in § 13.31 

3.L Breakout parameters 

The relation between the parameters at breakout and the 
physical parameters of the SN explosion were given in nu- 
merous publications ( e.g. M atzner & McKee 1999; Katz et al. 
|2009; Nakar & Sari 2010h . In particular, a complete set o f 
such relations is given in Appendix A of Nakar & S"^ (1201 Oh . 
For convenience, we reproduce the relation for vq and po be- 
low. As explained in § |2] these quantities completely define 
the planar problem. 

We use the approximate relation for the evolution of 
the shock velocity thr oughout the star (Eq. 19 in 
iMatzner & McKedl 19991) to set 

vo«1.0v4p/po]°■'^ (12) 

with 

v, = (£i„/Mej)'/2, p = Mej/(4^/?V3). (13) 

Here, Mej is the mass of the ejecta and Ein its energy (note 
that p is different from = M^jR-^ by a factor of 47r/3). 
We further use the density parametrization p{x) = fpj5(x/R)", 
where fg is a d imen sionless parameter of order unity (see 
ICalzavara & Matzneri 120041 appendix A for detailed esti- 
mates, note that pi/p* defined by the authors is related to fp 
defined here by fp = 47r/3(pi/p*)). Solving for r = /Jg', the 
following relations are obtained for n = 3 (appropriate for a 
blue supergiant (BSG)) and for n = 3/2 (appropriate for a red 
supergiant (RSG)): 

Vo/V* - 13A/io 8 5/^12 K0.4 fp [BSO) 

- zL '^il^013,,0.13 p-0.26 0.13 f-0.09 /nc/-N n A\ 
-4.5Mio 8 5/^13 K0.4 Jp (14) 

PO = 7 X lO-'KF^:%'lRiP'^^ffp''g cm-3 (BSG) 
= 2 X lO-'^MVh:%'!R-,l-'\^/'Y,'h cm-3 (RSG), (15) 

where Mej = lOMioM©, R = lO^^Rn cm = lO'^/^n cm, and 
V, = 3,000v,.8.5 km s"'. 

3.2. Applicability of the planar solution to a SN breakout 

Several effects must be taken into account when using the 
planar solution to describe the observed breakout burst. 

Light travel time smearing — A distant observer sees the break- 
out emission coming from different locations on the surface at 
different times due to the finite light travel time. This can be 
accounted for by appropriately "smearing" the instantaneous 
luminosity on a time scale fsmear ~ i?/c, as described in §|5] 



Relativistic corrections — At high velocities, /3o = Vq/c > 0.1, 
relativistic effects may introduce corrections of order tens of 
percents to the observed luminosity. One aspect which is eas- 
ily accounted for are first order corrections to the relation be- 
tween the calculated comoving fluxes and the observed fluxes. 
These corrections can be accounted for by Lorentz transform- 
ing the quantities in the comoving frame to the laboratory 
frame and finding the retarded time fret of the emission of pho- 
tons arriving at the observer at time fobs- This is discussed in 
§11 

3.3. Limitations 

The applicability of the planar solution is Umited by the 
following complications. 

Spherical expansion — As the outer mass elements expand, 
their optical depth decreases like r cx r~^, where r(m,t) is the 
radius to which the mass element moved. The planar approx- 
imation breaks once the optical depth of the outermost ele- 
ments drops significantly. Given that the outermost mass ele- 
ments move with a velocity v « 2v(), the optical depth of the 
outermost elements drops by a factor of ^ 2 at 

?sph^^^-^Wvo). (16) 
2 vo 

The use of the planar solution is limited to times f ^ fsph. 

For the planar solution to be applicable, is is required that 
fsph ^ to (this is equivalent to R^ xq). Using Eqs. (fT4l i-(fT6]l 
we have 

fo/fsph - QmM^,'\ti!R'l?-^."'C (BSG) 

^omM:r<,'^R'li^^^'f;"'' (RSG) (17) 

implying that fo <C fsph for practically all progenitors. 

We emphasize that our results are not applicable for pro- 
genitors with optically thick winds. 

Non spherically symmetric explosions — It is likely that SN ex- 
plosions are not spherically symmetric. The use of the planar 
solution for asymmetric explosions may be limited due to sev- 
eral effects: 

1. The break out velocities may be different at different 
locations on the surface; 

2. The shock may reach the surface at oblique angles; 

3. The shock arrival time to the surface may depend on 
location. 

The treatment of the first two effects is beyond the scope of 
this paper If the shock arrives to the surface at a large an- 
gle, the planar symmetry does not hold and our solution is 
not applicable. As the shock propagates through the star the 
anisotropy is expected to be smoothed, and it is reasonable 
to expect that there is a wide range of parameters for which 
the velocities are approximately uniform and the obliqueness 
is small. We note that even if the velocities are different, 
our solution can be applied locally, with the local value of 
the shock velocity, at any location where the shock arrives at 
small obliqueness. 

The third problem can, in principle, be considered within 
the context of the planar solution. It amounts to an appropriate 
smearing of the instantaneous luminosity over a time scale 
Afasym spanning the arrival of the shock to different locations. 
Section § |4] discusses properties that are not affected by this 
complication or by light travel smearing. 
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4. BREAKOUT ENERGY, MAXIMAL EJECTA 
VELOCITY AND ASYMPTOTIC LUMINOSITY 

In this section, robust approximate expressions are given for 
the total breakout energy, the velocity of the fastest moving el- 
ements, and the luminosity at late times, max(Afasyni,7?/c) < 
f < R/{Amq). These expressions are insensitive to light travel 
time averaging and to deviations from instantaneous arrival of 
the shock to the surface at all location. 

4. 1 . Breakout energy 

The breakout energy, defined as the total emitted energy up 
to the time of transition to spherical expansion, is given by 
(using Eq. (O) 



1- 



^sph 
fl,fO 



-1/3- 



(18) 



Using equation (fTTT l and a, ^ 0.1, the second term in the 
parenthesis is found to be of order 0.1 and to a good approxi- 
mation ^bo ~ £00 ■ 

The breakout energy is not sensitive to the precise value 
taken for the time of transition to spherical expansion, as we 
now illustrate. Beyond the transition to spherical expansion, 
the luminosity will level off approaching an asymptotic power 
law 

L(f)aisph(fAsphr"^''% (19) 



the luminos ity at f^ ph and 
7,(312) (iChevaheii IT992I: 



where Lsph = L(/sph) is 

a, pi, = 0.34(0.17) for n 

iRabinak & WaxmanI I20T0I: iNakar&Saril 120101) . It is 
useful to estimate the time at which the contribution of the 
later spherical phase emission, A/s, is comparable to E^o- 
The time it takes to accumulate an energy E^o in the spherical 
phase is roughly given by 



AfBO = fsph [(1 -asph)£'Bo/(isphfsph)] 
/ t \ l/(3-3asph) 

\a,tQ J 

Using Eqs. (O, ([17]) and ([T4| we find 



l/(l-a.ph) 



(20) 



AfBov,//?= 1.8M■ 



-o.ol,,-o.ol d003 



10 



Rn'^^'^A'^fp'' (BSG) 



= 2.4M?n''^?f.«7?'4r 



i.5 

'10 8. 5^13 

implying that for all progenitors 



f/' (RSG), (21) 



(22) 



At earlier times, the accumulated energy emitted in the spher- 
ical phase, AE, relative to the breakout energy, is roughly 
given by 

— - (^-] (23) 

EbO \ AfBO / 



The difference in shock arrival time to the surface due to 
asymmetry is always much shorter than /?/v*, implying that 
the breakout energy can be accurately integrated, even if the 
time difference in arrival times is significant. 

4.2. Maximal velocity 

The velocity of the fastest moving ejecta can be obtained 
from the planar relation Eq. (fTOt . This velocity can be probed 
by other observations, including the spectrum of the breakout 



and radio observations of the collisionl ess shock propagating 
through the circumstellar m edium (e.g. lChevaher & FranssonI 
I2OO6I: IWaxman et"an 12007 ). As long as the planar approxi- 
mation is vahd, the velocity of the surface is proportional to 
the emitted energy and is given by Eq. (fTOl i. Once the ra- 
dius changes considerably, the acceleration declines sharply, 
d,v oc LR~^, and th e velocity does not increa se significantly 
any more (see also ' Matzner & McKeelll999l) . The resulting 
velocity of the fastest part of the ejecta is thus related to the 
breakout energy 



47ri?2, 



-Ebo, 



and is given by 



:2.0V0 



^sph 
fl,f0 



-1/3' 



(24) 



(25) 



A rough estimate of the additional acceleration beyond the 
transition to spherical expansion can be obtained by approxi- 
mating L(t > fsph) = const = Lsph and R oc t. In this limit, the 
acceleration drops like f beyond fsph and the total additional 
velocity Av^ax is roughly given by [using Eq. ( fTsT ll 



Av„ 



^sph^sph 

Ebo 



fsph 

a, to 



-1/3 



(26) 



implying a correction of a few percent for all progenitors con- 
sidered. 



4.3. Asymptotic luminosity 

In the regime fsmear < t < fsph, where fsn,ear = 

max(Afasym,^/c), the luminosity is given by equation 
([8]l and can be expressed as 



L(f) = ioo(fAo)"^/^ 
-nR^K 



4 -2-4/3(:^ 



3.0 



1/3 



,.4/3.-4/3 



= 2AxlO''R],n-^r<9Po/k:''^^ES-\ 



(27) 



where vo.g = lO^vo cm s"', po,-9 = lO'^pog cm"-' and f = Ifhr hr 
The weak dependence on the parameters po and vq implies 
that, if detected, this power law tail can be used for an accurate 
determination of the stellar radius. 

We emphasize that since f smear >R/c and fsph ^ R/4-vo, the 
time interval fsmear ^ f ^ fsph exists only for small velocities 
satisfying vo/c <C 1 /4. This condition is met only for large 
RSG progenitors. Given that vp /pp oc R^'^^ for n = 3/2, it is 
useful to rewrite the relation (|27) as 



Rv. 



4-71 I 



Vo 



1 o inI3r2/5 S/15 

= 1.8 X 10 L43 fj^. cm 

/ 3/2\ ^/'-^ 

\ Vo,9 / 



(28) 



where L = 10^^L43 erg s ' and the factor in the last line is close 
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to unity for RSG parameters. Using Eqs (fl4t and (ITSt . 

/ „3/2\2/15 



V0,9 



. 1 1 ^0.{)26„-{).24p0.{)15 ^0.43 ..0.07 (pcr\ 



(29) 



5. LIGHT CURVE 



In this section, the observed bolometric light curve is cal- 
culated assuming that the breakout is strictly spherically sym- 
metric (i.e. assuming negligible difference in shock arrival 
times to the surface at different locations). In § 15. II e xact ex- 
pressions for the observed luminosity are given. In § 15.21 an 
approximate analytic expression (Eq. d3&t ) is derived assum- 
ing cti)/R ^ 1, which is valid for all progenitors with the ex- 
ception of the largest RSGs. Approximate expressions for the 
peak luminosity and the luminosity at the time R/c are given. 
The results of this section are summarized in § 15.31 

5.1. Exact light curve 

Due to the difference in arrival times of photons originating 
from different positions on the surface of the star, the actual 
luminosity Lobs(f), that a distant observer would measure, is 
related to the instantaneous luminosity L(t) by 



■^'obs(f)= / h(n)L(t-R(l- n)/c)^id^i, 



(30) 



where h(ii) = 2TTl(p)\r=Q/ C is the angular distribution of 
the radiation intensity at the surface normalized so that 

^|^h{p)^ld^l=\. 

The precise value of h(p) requires the solution of radiation 
transport up to the surface. For the non-relativistic breakouts 
considered here, the transport equations can be solved in the 
steady state approximation with the flux given by the diffusion 
solution. For the considered c ase of Thomson scattering, h{pL) 
was obtained analytically by IChandrasekhad (I1950D and can 
be fit by a linear relation. 



with 



h{^) ai + bj^, fl//2 + /7//3 = 1, 



fl/ = 0.85, /7/= 1.725, 



(31) 



(32) 



to a good approximation (better than 3%, see §0. For com- 
parison, Eq. dSTT i with aj = 2, bj = represents a black body 
surface (isotropic emission), while a/ = 1 , hi = 3/2 represents 
isotropic scattering in the Eddington approximation. The re- 
sulting light curves for a density power law index n = 3 and 
for different values of cto/R are plotted in figures[T]and|2] 

5.2. Approximate light curve for R/ c 3> ?o 

For most progenitors, the smearing time scale R/c is much 
larger than fo. In fact, using Eqs. ( fT4b and ( fTSl l we have 

ct,/R = Q.Q2MT'vt^!R^!n,TC (BSG) 

= 0.05Mrr'v;'8'|^l3'^a4''./;"'' (RSG). (33) 

As can be seen, except for very large RSG's with R ^ lO'"^ cm, 
cto/R may be assumed to be small. 

In this case, the burst time scale is R/c and the typical lu- 
minosity is Eoqc/R. It is useful to describe Lobs as a function 
of 

C(t — f|-ef ) 

s=^-^. (34) 



Consider first the formal limit t{)c/R —?' 0. In this limit L(t) — 
EooS{t-ti-gf) and the observed luminosity ( |30] | goes to 



Lobs(0 



tac/R^i) R 



■{\-s)hi\-s)-{0<s<\). (35) 



In figure [T] the light curves are shown as a function of s 
(blue solid lines) while the limit of Eq. (|35] | is shown for com- 
parison (dashed red). As can be seen, the observed luminosity 
converges slowly with cfo /R to the limiting value. This is due 
to the luminosity tale of L oc f~^/^ at late times, the integral 
of which converges slowly. This can be taken into account by 
using the approximation of Eq. (|8]l. Using Eqs. OTl i. (l30l l and 
^ we find 



Lobs(0 = 



where 



CEr 



R 



1 pmin{\ .s-Sc) 



3 7o 



(\-s')[ai+bi(\-s')]{s-s')-^l^ds' , 
(36) 

Sc = catto/R. (37) 



An e xplicit expression for the integral in (l3Ft is given in Eq. 
( IA2I ) and the resulting observed luminosities are shown in fig- 
ure |2] (dashed dotted magenta lines, values of Eoo/Eq = 2.03 
and fl, = 0. 1 fitted for the case n = 3 considered were used). As 
can be seen in the figure, this is an excellent approximation to 
the calculated observed luminosity. 

A few properties of the light curve can be derived from Eq. 
( [36] l. At very large 5^1 Eq. ( [36] l reduces to equation ( |27T i 
as required. Given that ic ^ 1, the value of Lobs at i = 1 (f = 

tret + R/c) is 



Ubs{t = tref + R/c)-- 



cEc 



R 



s'JHa,/2 + b,/5), (38) 



where for a, = 0.85(2) we have (a,/2+b,/5) = 0.77(1). Com- 
paring to ( |27] i (and using the relation between a,,Eoo,Loo in 
(|9]l), we see that at t = R/c the luminosity drops to a value 
about 3 times larger than the extrapolation of the asymptotic 
luminosity (Eq. ( l27b ) to this time. 

As shown in §|A] the peak observed luminosity is, to a good 
approximation, given by (Eq. (|A5t ) 



Lpeak = (a/ + ^/)- 



R 



1- 



1/4- 



(39) 



5.3. Light curve calculation summary 

The observed light curve can be calculated usi ng Eqs (l30i 
and (ISTT i. with L(f ) given by (|6]l and tabulated in Sapir et aLl 
( 1201 lal . table 3). The following parameters completely deter- 
mine the light curve: breakout luminosity Lq, breakout shock 
crossing time fo, time of peak planar flux f,ef , proge nitor radius 
R and d ensity power law index n. As shown in ' Saprr et al.l 
(1201 lal) . the luminosity depends weakly on n. The parame- 
ters Lq and fo are related to the breakout density, velocity and 
opacity through Eqs. (|7]i and (|5]l. The resulting light curves 
for the case n = 3 are shown in figure |2] 

For most progenitors R/c^ t^ and the observed light curve 
can be calculated using Eq. ( |36] | with s given by Eq. ( l34b and 
ai = 0.85, bi = 1.725 given by E q. ( l32b . An explicit algebraic 
expression is given in Eq. (IA2b . In this case the light curve 
is determined by the following parameters: breakout energy 
Loo, time of peak planar flux fjet, progenitor radius R and a 
dimensionless parameter Sc- The parameters Loo and Sc are 
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Fig. 1 . — Normalized observed luminosities as a function of the parameter 
i = c(f — ;,-ef)/S for different stellar radii (blue solid lines), compared to the 
limit of Eq. (35) (dashed red). 




Fig. 2. — Normalized observed luminosities (solid blue lines) as a function 
of time since (expected) breakout, for different star radii . Also plotted are 
the instantaneous emitted luminosity (bl ack dashed line) and the analytical 
observed luminosity given by Eqs. l |36) and | |A21 (magenta dashed-dotted 
line). 

related to the breakout density, velocity and opacity (and R) 
through Eqs. (01 and 

The peak luminosity in this case can be approximated by 
Eq. ( [39] l. At f = tref + R/c the luminosity drops by a factor 
^ (O.lcto/Ry^^ compared to the peak, to a value which is 
approximately 3 times higher than the extrapolation (to t = 
trei+R/c) of the asymptotic luminosity given by Eq. ( |27] |. 

In the extreme limit fo — > a simple approximation for the 
light curve (for t < R/c) is given by (l35T l. which depends on 
two parameters only, £00 and R. 

6. COMPARISON TO PREVIOUS STUDIES 



In § 16.11 we compa re our results with those of the nu- 
merical calculations of lEnsma"n &Burrowsl([T992l) for 1987 A 
like (BSG) progenitors. In § I6.2l we compare our calibrated 
analytic results w it h the order of magnitude estimates of 
iMatzner & McKeel ( Il99l and lNakar & Saril ilOld) . 

6.1. Comparison to Wnsman 



In figu re[3]we compare our results with the numerical light 
curves of lEnsman & Burrowsl (11992b . As can be seen, the nu- 
merical light curves are in excellent agreement with our an- 
alytic light curves, for an appropriate choice of po and vq. 
The fitted values are {vq = 16,500 km s"' , po = 1 .6 x lO"''} 
and {vo = 25, 000 km s~' , po = \.2x lO"**} for the progenitor 
models '500fulir and '500full2' respectively (the opacity is 
assumed to be k = 0.33 cm^g"', appropriate to a mixture of 
ionized Hydrogen and Helium with mass fractions X = 0.67 
a nd Y = 0.33). 

lEnsman & Burrowsl (119921) have used e + p = 2.5C/c for 
the calculation of the instantaneous lab frame luminosity, 
Lvdh = 4-nR^[C + v(e + p)], which is not accurate to first order 
in /3 (see §|B]for a careful calculation of the frame transfor- 
mation), and assumed isotropic emission ih(fi) = 2) for the 
light travel time smearing calculation, which does not repre- 
sent Thompson scattering opacity (see § 15. 11 1. For the com- 
parison we have adopted similar values in our calculation of 
the light curves shown in figure |3] although they do not yield 
an accurate description of the light curves. 

Using Eq. ( fT6b . the transition to the spherical phase is ex- 
pected at approximately 7400 s an d 7100 s respectively (usin g 
the reference time as in figure 2 in l Ensman & Burrowslll992"l) . 
As can be seen in figure [3] the planar solution agrees with the 
spherical calculation at the latest times presented (which are 
less than fsph) up to a factor of ^ 1.5. The maximal veloc- 
ities predicted by equation Eq. ( |25] l are 28 , 000 km s~' and 
43,000 km s"', in good agreement (up to ^ 10%) with the 
maximal velocities of 30,000 km s"' and 48,000 km s~' ob- 
tained in lEnsman & Burrowsl (i 1992b . 



45 



; 2.3(51) 



43 



3. 
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6700 S800 6900 7000 7100 7200 
time for 10^' erg mode! (sec) 



7300 



Fig. 3. — A compari son of our light curves with those obtained by the nu- 
merical calculations of lEnsman & BurrovifsH1992fl for two progenitor models 
('50()fulir and '500full2'). Our results are shown in color (blue and red for 
source frame a nd observed luminosities r espectively), overlaid on the oiigi- 
nal figure 2 of lEnsman & Burrowsl dl992D. For the comparison we have used 
e+p = 2.5Z^/c and hi/j.) = 2, as used bv lEnsman & Burrows ( 1992), although 
these values do not yield an accurate desciiption of the light curves (see text). 



6.2. Comparison with previous analytical work 

A quantitative comparison of the analytic result s presented 
here ( calibrated using the numerical results of ISapir et all 
1201 lab to those of the order of magnitude estimates made in 
the literature is compUcated since the breakout parameters po 
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and vo were not well defined in earlier analyses. We remind 
the reader that we define po and vq to be the upstream density 
and the shock velocity that would be obtained in a pure hydro- 
dynamic calculation (ignoring radiation diffusion) at the point 
where t = c/vj (see §|2|. For illustration purposes we com- 
pare a fe w of our results with orde r of mag nitude estimates 
made by jMatzner & McKedl 19991) and ( N akar & Sarill2010l) 
in which po and vq were (vaguely) defined in a similar way. 

• The numerical values of the order of magnitude es- 
timates of the power law decUne of the luminosity, 
L cx f"^/^, made by dNakar & Sarill2010l) can be summa- 
rized as L(t) - 47r7?2K-4/3(vo/po)'/^c4/3r'*/3 (using the 
values for po and vq given in their appendix). This re- 
sult is in rough agreement with our Eq. (|27| . with their 
expressions overestimating the luminosity by a factor 
of about 3. 

• iMatzner & McKee' (1999) estimate that the maximal 
velocity obtained by the fastest ejecta is about twice the 
breakout shock velocity vq. Coincidently (as far as we 
can tell), this turns out to be accurate (see §l25t. 

• The numerical values for the breakout energy of RSGs 
and BSGs in ( Matzner & McKee 1999) are related to 
the numerical values of the maximal velocities by Ebo = 
47r/?^cVmax/K and Ebo = 1.7 x ATrR^cw^nax/ k respec- 
tively. The RSG relation is accurate while in the BSG 
relation the emitted energy is overestimated by a (mod- 
est) factor of about 1 .7 (see (l24b). 

7. SUMMARY & DISCUSSION 

We have derived exact bolometric light curves of supernova 
shock breakouts using the universal planar breakout solu- 
tions (Sapiret al. 2011a), assuming spherical symmetry, con- 
stant Thomson scattering opacity, k = (Z/A) aj /nip, and an- 
gular intensity corresponding to the constant flux limit. The 
light curves are insensitive to the form of the density profile. 
This was demonstrated by calculating the emission for power 
law profiles p cx x" with resulting luminosities changing by 
< 30% for a broad range of power law indexes, 1 < n < 10 
(ISapiretal.l | 2011al) . 

The breakout emission properties are determined by four 
dimensional parameters: the progenitor radius R, the break- 
out velocity and density, vq and po respectively, and n. vq 
and pq are the shock velocity and (pre-shock) density at the 
point where = c/vsh is reached in the pure hydrodynamic 
(neglecting radiation diffusion) solution (see § |2] for exact 
definitions). The relations between the SN parameters, the 
ejecta mass Mej and bulk velocity v,, and the breakout pa- 
ram eters, pn and vp, a re given in Eqs. (fT4l i and (flST l (see 
also IMatzner & McKed [T999: Katz et al. 2009; Nakar & Sa3 

Hoim 

The application of the planar solution to SN breakouts 
was discussed in § |3] The planar approximation is applica- 
ble provided that the shock crossing time at breakout, fo = 
c/kpoVq, is much smaller than the time for tr ansition to spher- 
ical expansion, fsnh ^ ^/(4vo) (Eq. (O, Piro et aLlfioiOl: 
[Nakar & Sari 2010). This is valid for practically all progeni- 
tors, see Eq. (fTTI l. At f > fjph the expansion is no longer planar, 
and the planar approximation no longer holds. Our results are 
not applicable for progenitors with optically thick winds. 

Deviation from spherical symmetry may affect breakout 
light curves due to several effects: the breakout velocity may 



be different at different surface locations, the shock my reach 
the surface at oblique angles, and the shock arrival time at the 
surface may depend on location. In this paper we have as- 
sumed that the shock reaches the surface with approximately 
the same velocity everywhere and parallel to the surface. Dif- 
ferences in shock arrival times to the surface, Afasym, due 
to moderat ely asymmetric explosions, affect the early light 
curve (e.g. Calzavara & Matznerll2004b . at f < Afasym, but do 
not affect the maximal velocity of the ejecta and the total emit- 
ted energy, v^ax and Ebo- 

Analytic expressions for Vmax, Ebo and the late time lumi- 
nosity were derived in section § |4] The total energy of the 
breakout burst is approximately (Eq. (fT¥t ) 

Ebo = SttR^k-^vqc = 1 .9 x lO"*^ erg^^jVg^o.V (40) 

where vq = lO^vg cm sec"'. We have shown that integrating 
the luminosity to times greater than fsph affects the total energy 
considerably only at very late times, t >R/v^, (Eq. ^ZTV ). The 
maximal velocity of the ejecta is directly related to the emitted 
energy by (Eq. (|24] |. see also lLasher & Chanll979HSapir et al.l 
l20Tlal) 

Vmax = KEBo/(4TrR^c) = 2 .Ovq . (41) 

We note that these results are not sensitive to deviations from 
of the steady state flux angular intensity distribution we have 
used. For 4vo ^ c, valid for large RSG progenitors, there 
is a significant separation between R/c and fsph- In this case, 
the luminosity at max(Afasym,^/c) <t < fsph is approximately 
given by (Eq. ( |27] |) 

L(f) = (4/3)7r/;2(vo/p„)'/3(c/K)4/3f-4/3. (42) 

The strong dependence of the asymptotic luminosity on R and 
weak dependence on po and vq allows one to accurately deter- 
mine the progenitor radius of RSG breakouts (Eq. (l2Ft ). 

R = 2xlO"L%^tll^\m. (43) 

The bolometric light curve, assuming negligible spread in 
shock arrival times, are calculated in §|5] A proper calculation 
of the effects of finite light travel time requires knowledge 
of the angular dependence of the intensity. Fortunately, the 
problem of radiation transport in an optically thick medium 
with opacity do minated by Thomps on scattering was solved 
in closed form ( ChandrasekharllT950 . results summarized in 
§ 0. Exact light curves can be calculated using Eqs. ( |30] | 
and {T\\ and the planar luminosity functions, C, given in 
(lSapiretal.ll2"0TTa h. Some examples are shown in figures [T] 
and|2] For cases where cfo ^ R, applicable in all progenitors 
except for the largest RSGs (see Eq. l33T l. the planar luminos- 
ity C can be approximated by a power law, Eq. (|8), allowing 
the derivation of an analytical expression for the light curve, 
given in Eqs. (|36] | and (IA2b . The analytic expression is com- 
pared to the exact calculation (both without relativistic cor- 
rections) in figures [T] and |2] In this case, the peak luminosity 
is approximately given by 

iobs,peak = 2.5iEBoc/R)ll-icto/R)'^\ (44) 

where 

Eboc/R = 5.6 X 10'''^7;i3V9Ko'4 ergsec"' (45) 

is the typical peak luminosity. In addition, we have shown that 
at t = R/c the normalized luminosity Lf"*/^ is about 3 times 
larger than its asymptotic value given by Eq. (l27t . A short 
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summary which explains how to use the diffe rent expressions 
to obtain the light curves is provided in § 15.31 

The transformation of the rest frame intensity to the lab 
frame and the value of the retarded time introduces correc- 
tions of order /?'. These are calculated in §|B] For the veloc- 
ities considered, /3o ^ 0.3, the only considerable correction is 
to the early light curve and is given by ( IB 19) . We note that 
other corrections of order /?' are not excluded. 

The results of this paper are compared to previous results 
in § ID The calculated bolometric light curve is shown to 
be in excellent agreement with the numerical calculation of 
lEnsman & Burrowsl (119921). see figure E The order of mag- 
nitude estima t es giv en by iMatzner & McKed (119991) and by 
iNakar & Saril (120101) for the emitted energy, maximal velocity 
and asymptotic luminosity agree to within factors of few with 
our exact analytic expressions. 

The properties of the breakout emission depend strongly on 
the radius of the progenitor R and on the breakout shock ve- 
locity Vo, depend weakly on the value of the density at break- 
out pq, and are insensitive to the density structure. Break- 



out observations therefore allow one to accurately determine 
R and vq. These quantities are directly related to other ob- 
servables: The breakout shock velocity vq is roughly pro- 
portional to the ejecta velocity v,, which is probed by SN 
observations. The maximal velocity of the ejecta (about 
twice the breakout shock velocity) can be constrained by 
the radio and X-ray emission produced by th e interaction of 
the ej e cta with the circumstel lar medium (e.g. IWaxman et al.l 
I2007t ISoderberg et"ari 120081) . Both parameters affect the 
subsequent spherical expansion (cooling envelope) phase of 
the emission (e.g. C hevalier.. 1992: .Rabinak & Waxmaa.201C)l: 
INakar & Sarill2010h . 
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APPENDIX 
A. LIGHT TRAVEL TIME AVERAGING 



The integral 



1/3 i 



min(l ,^-5^-) 



(l-s')[a, + bi(l-s')](s- s'r^l^ds' 



(Al) 



appearing in equation ( [36] l is explicitly given by 

{l^y ' [-ai{s - 1 )2/3 +6-bi{s-\fl^-[ai{\-s) + 2bi( i - s + fs^ )] j-i /3] s>\+s, 
a,{\-s-\sc) + bi [\ + (s,-2)s + s^-Sci\ + ^Sc)\ +isc/s)^l^ [-a,-b,+ l{a, + 2bi)s-lbis^] <s<l+s, (A2) 

S<Sc 

To lowest order in s^, the peak of g(s;sc) is reached at 



with a peak value of 



gpeak= (fl/ + ^/) 



l_(3-3/4 + 3l/4)('^£l^ 

V ai + b; 



1/4 



(A3) 



(A4) 



Using Eq. the fact that for < a,<2 we have 1 < [(a,+2b,)/(a,+b,)]^/'* < 1.2, and the numerical value (3"^/"^+ 3 '/'*)a''''* = 
1.0 (using a, = 0.125, see Eq. (|9]l), we conclude that to a good approximation 



iobs.peak - — —ifll + bi) 
R 



1/4' 



(A5) 



B. FIRST ORDER CORRECTIONS DUE TO THE TRANSFORMATION BETWEEN THE COMOVING FRAME 

AND THE OBSERVER FRAME 

In this section the first order corrections in /3 due to the transformation between the comoving frame and the observer frame 
(lab frame) are calculated. Note that the expression for the asymptotic luminosity ( |27l ) is valid only for small velocities vq ^ c/4, 
for which the first ord er co rrections are negligible. We estimate the first o rder corrections in /3o = Vo/c to the value of the breakout 
energy, Eq. (fTsT l, in § IB.1I and the corrections to the light curve in § IB.2I 

Bl. First order corrections to the Breakout Energy 
The rate of energy escaping from the surface of the expanding envelope as measured in the laboratory frame, dE\^\, / dt, is given 



by 



dEu 



dt 



= 47rr (Aab-veiab), 



(Bl) 
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where r(f) is the radius of the surface and £iab and eiab are the radiation flux and energy density at the surface in the laboratory 
frame. The values of £iab and eiab can be expressed using the comoving frame flux C, energy density e and pressure p by 



£iab = /: + v(;? + e) + (9(/32), 

eiab = e + 2/3c-'/: + 0(/32). (B2) 



The emitted energy per unit time to first order in /3 is thus 



T=() 



C. (B3) 



The term {pc/C)r=o depends on the angular distribution of the radiation intensity at the surface, which dep ends on the transport 
properties of the medium. For Thomson scattering in the constant flux limit, the value is (pc/C)r=o = 0.71 tChandrasekhaJ 1 93oi 
see also §|C]'. For comparison, for a black body surface (isotropic emission), the value is (pc/£)r=o = 2/3 and for isotropic 
scattering in the Eddington approximation it is {pc/ C)r={) = 17/24. 

Assuming that {pc/£)t=o is constant over time and neglecting the change in the surface radius during the emission, Eq. (IB3t 
can be analytically integrated over time. To do this note that AttR^ J C{t')dt' = E(t) and that 

/' c f' dE 1 

C(t')v{t')dt' = - — Edt' = -E(t)v(t), (B4) 
Kjdt' 2 

where Eq. ( fTOl l was used. Using Eqs. (IB4b the integration of equation iB3i to infinity results in 

^lab.oo = £00 f 1 + X ^ Poo], (B5) 
\ 2 L T=0 J 

wher e Po p = Voo/c is the asymptotic value of /3 in the planar approximation (approximately unchanged by spherical geometry, 
see § l4.3l l. Using the approximation of Eq. (IB41 i. and adopting {pc/C)t=o = 0.7, we have 

(l + 0.7/3o) = £oo(l+0.35/3oo). (B6) 

Linearly adding this correction and the small term (hph/(a,to))~^^^ in equation ( fTSl l we obtain 

-1/3 



l-(^) +0.7/3o 



(B7) 



B2. First order correction to the observed Light curve 

First order corrections to the light curve arise from corrections to the intensity, direction and value of the retarded time. At 
times t > R/c, the surface of the star moved a distance of {m /c)R implying that there are corrections of order v/c due to the 
spherical nature of the expansion. An estimate of this correction is beyond the scope of this paper Here we focus on the first 
order corrections at early times t R/c (including the peak observed luminosity). At these early times the angle between the 
emitting region radius and the direction towards the observer is small, /i- 1 <C 1, allowing a simple derivation of the first order 
correction. 

Consider a spherically symmetric moving surface with radius r(t) emitting radiation with lab frame intensity I{t;il) which is 
axisymmetric with respect to the surface normal. Consider an observer located in the direction fiobs at a very large distance. The 
position on the sphere is parametrized by the angle 6 between il and fiobs- Since the intensity / is assumed to be the same at any 
position on the sphere and to be axisymmetric with respect to the normal we have /(fjflobs) = l{t,0). The luminosity inferred by 
the observer is given by 



iobs(fobs) = 47r y 2TTadaI(t;e), (B8) 

where 

a = rsin6' (B9) 
and the lab (retarded) time t is related to the observer time fobs through 

R + ct-rcos9 = ctobs- (BIO) 

The integral (|B8t should be evaluated at constant fobs, with f and 9 functions of a through the relations ( |B9t and (IB 1 Oi l, in the 
regime (3 < cosd < 1 (it is assumed that there are no photons with cos 6 < P, which are coming from outside of the surface). 
It is useful to solve for /i = cos 9 

i?-c(fobs-f) 

p= , (Bll) 

r 

where f' = fobs-?- Using (IB 10b . we find 

cdt = rdji + jjNdt => rdji = c( 1 - (iij)dt (B 12) 
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and 

ada = r sin^ 9vdt + sin 6 cos 9d9 

= cr[il-i/)P-Hil- Pii)]dt = criP-n)dt. (B 13) 

Equation dBSI I can be written as 

i'obs(fobs) = 47r J dtcrip-P)2^TI(t■,^i) (B14) 

with /i given by dBl II ). 
The lab frame fi is related to the surface frame /i, /icom, by 

Aicom=-; — (B15) 
while the lab frame intensity / is related to the surface frame /com by 

/(f , /^) = [7( 1 + /?Mcom)] ^/com(f , /^com) , (B 1 6) 

where 7 = (1 - /3^)'/^ is the surface Lorenz factor 

Assuming that the angular dependence of the comoving intensity is time independent, /com (f, A*) = KlJ)je,com(t)/('2Tr), equation 
( IB 14b can be written as 

Lobs(tobs) = J dt^(p-l3)'y^(l+/3^i^om?h(pcom)Lcom(t), (B17) 

where fi and fi^om are given in eqs (IBl lb and ( |B15b . 

Consider next the first order corrections in the parameters /? and (1 neglecting the difference between r and R. Note that 
in this approximation we can set /J/x — >■ /5, 7 —> 1 , /icom — > M ^^d 

iobs(fob.s) = y" dtfih(p)( 1 + 2/3)Lcom (B 1 8) 

with^= l-c(fobs-f)/^- 

Using L{t)\{t)dt « jVoo/ioo ~ 2vo£'oo (see (IB4I )). we find that in the limit t()c/R ^ 1, the correction at early times fobs -^R/c 

is 

iobs=iobs,o(l+2A)). (B19) 

C. STEADY STATE RADIATION TRANSPORT WITH THOMSON SCATTERING 

The steady state pr oblem of radiation transfer in a semi-infinite medium with Thompson scattering was analytically solved by 
(IChandrasekharl[T950) . Assuming a constant flux C coming from inside the medium and zero incident radiation, the intensities 
at the surface, /;(/j.) and /, (//), polarized in and perpendicular to the meridian plane respectively, were found in closed form. The 
polarized intensities are given by 

^' = ;|^^^'^^^' 
1 3 

I,-= ^—CHr(p)(p + c), 

(CI) 

where //, (; = l,r) are the solutions to the integral equations 

//,(/i) = 1 + f^Hii^i) [ ^i^lzlp.H,(^')dfi' (C2) 
Jo M + M 

with a; = 3/4 and = 3/8, while q and c are the solutions to the equations 

q^ = 2(l-c^), 
qHiil) = (l+c)HXl). 

(C3) 

The resulting intensities were numerically calculated by Chandrasekhar and his secretary, Mrs. Frances H. Breen, to 5 decimal 
points for the twenty values of /i = : 0.05 : 1 using pen and paper (Table XXIV in .Chandrasekhar. 1950. note that F = C/tt). The 
resulting total intensity, / = // + /r, can be fit by the Unear relation 

hch{^i) = ^ « 0-85 + 1 .725/^ (C4) 
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to an accuracy better than 3% for all < /i < 1 (by assumption, /(/i) = for /i < 0) . 
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